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Abstract 

Emergent anyons are the key elements of the topological quantum computation and topological 
quantum memory. We study a two-component fermion model with conventional two-body inter- 
action in an open boundary condition and show that several subsets in the low-lying excitations 
obey the same fusion rules as those of the toric code model. Those string-like non-local excitations 
in a given subset obey mutual semionic statistics. We show how to peel off one of such subset 
from other degenerate subsets and manipulate anyons in cold dipolar Fermi atoms or cold dipolar 
fermionic heteronuclear molecules in optical lattices by means of the established techniques. 
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Introductions Anyons are long wistful objects in two-dimensional condensed matter systems 
. Interesting to explore anyons is currently renewed because of the potential application 



of anyons in topological quantum computation and topological quantum memory 
researches mainly focus on two topics: non-abelian fractional quantum Hall states 5|, [6 
and particular lattice spin models, e.g., Kitaev toric code model [3j, Levin- Wen model 
and Kitaev honeycomb-lattice spin model j^]. 
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Experimentally exciting, manipu 
or tried for the toric code model 10 



ating and detecting abelian anyons have been suggested 



ll| and for the insulating phase of Kitaev honeycomb- 



lattice model [121 ] . Although these non-trivial tries are interesting, the reliable evidence 
for the existence of anyons still lacks since neither the toric code model nor honeycomb- 
lattice spin model is easy to be realized owing to those unconventional interactions between 
constitution particles, e.g., cold atoms [13]. 



A widely interesting question is: do we have a simple system with conventional two-body 
interactions whose non-local elementary excitations are anyons? We here study a two- 
component fermion system which is a set of decoupled Ising chains 14[] . These chains form a 
two-dimensional square lattice with each chain along the horizontal diagonal direction (See 
Fig. [p. The fermions within a chain interact with on-site repulsive and nearest neighbor 
attractive potentials. Nussinov and Ortiz [15| have found the decoupled Ising chains are of 
the same spectrum as the toric code model and Liven- Wen model . However, the topological 
order is trivial when the plain lattice is curled up to a torus. 

We here show that although the topological order of decoupled Ising chains with a peri- 
odic boundary condition is trivial, when the in-chain coupling constants are special chosen 
with respect to the chemical potential, there are low-lying non-local excitations in an open 
boundary condition which may obey anyonic statistics. We will see that the low-lying ex- 
citations of this model are classified by several kinds of closed subsets. One kind of them 
is local, including a single hole, a double occupant and a spin-flip. Other two kinds consist 
of a local fermion and two string-like non-local excitations. The fusion rules and exchange 
phase factor in later two are exactly the same as those of the excitations in the toric code 
model. The statistics of the non-local excitations themselves are bosonic while it is mutual 
semionic. We find that the ground state in this model is stable for a dipole-dipole long range 
interactions if inter-chain couplings are negligible. Thus, this system can be realized in cold 
dipolar Fermi atoms, e.g., rare-earth atoms of Ytterbium 161 ]. or the cold fermionic het- 
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eronuclear molecules like 40 K 87 Rb [17], in optical lattice. With respect to the ground state 
degeneracy, any subset of the excitations is accompanied by many energetically degenerate 
subsets. We discuss how to peel off a subset from other degenerate subsets and manipulate 
anyons by means of the established cold atom and molecular techniques. 
Two-component fermion model in square lattice We consider a simple Hamiltonian for two- 
component fermions in a square lattice (Fig. [I]) 

H = - J.(2n Bii -l)(2n ad -l) + U y %2(2n hi -i)(2n lii -i), (1) 

(ij)hd,s i 

where n St i = c\ ^ S) , with c S) j being annihilation operators of (pseudo)spin-s fermions. The 
symbol (ij)hd means the sum is over nearest neighbors along the horizontal diagonals of 
squares. This is a set of decoupled Ising chains along the horizontal diagonals of squares. 
We restrict on the nearest neighbor attractive interaction while on-site is repulsive, i.e., 
J s > and U > 0. In the case, the ground states of this Hamiltonian are 2 n -fold degenerate, 
i.e., every individual chain is ferromagnetic, i.e., for a set of spins {s\, ■ ■ ■ , s a , ■ ■ ■ , s n }, 

n 

iG , w)=ni G -)=n c Ui°) ) (2) 

a=l a,2 a 

where n is the number of chains, i a is the site index in the a-th chain. Restricted to the 
open boundary condition, the low-lying excitations above a given ground state | for a 
given site i a , reads (See Fig. 1 for a given ground state \G^) which will be defined later.) 

Ki a \G {s} ) = (cj a>ia + c Satia )\G {s} ) = c Satia \G {s} ), 

AJG {S} ) = (cl M + c- Sa ,aG {s} ) = cljG {s} ), 

Fia\G{ s} ) = m ia V ia \G {s} ) =ic Sajio 4 a) JG {s} ), 

Wp\G {s] ) = U^ b \G {s} )= n</> c L'J G M>, 

W H p,p>\G {s} ) = J] ^H ta \G {s} ) = J] ic^^c^lGw), (3) 

w d P , P <\G {s} ) = n^*-i G w>= n ^A,>4 a jG {s} ), 

where P and P' denote two plaquettes on the right and left of i a , respectively, s =| (|) 
if s =| (|). The order of sites is defined by i' b < i a if i' b is on the left hand of i a for 
b — a or i' b is in a chain lower than the chain with i a . TljVjJ 7 create a hole, a double 
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occupant, and a spin-flip. W, W d and W h create a half-infinite string of spin-flips, a spin- 
flip string with a double occupant and a spin-flip string with a hole, respectively, since the 
spins of fermions at sites i' b < i a are flipped from their ground state configuration while 
those at j c > i a keep in their ground state configuration. The excitation energies of these 
local and non-local excitations in turn are 4J Sa + 2U, 4J 5a + 2U, 4Jj + 4Jj_, 2Jj + 2J|, 
2J| + 2Jj_ + 2U and 2 Jf + 2J^ + 2Z7", respectively. The finite energies of the string-like 
excitations mean they are deconfined. Note that the ground state may be instable against 
an on-site energy V B Y^i(^ n t(l)4 — 1) if > an d if Vp < 0, a pair of string-like excitations 
may be linearly confined and the statistics of the string-like excitations becomes ill-defined. 
These excitations are highly degenerate due to the degeneracy of the ground states. For 
example, if {s 1 , ■ • ■ , s a ■ • ■ , s n } -> {s 1 , ■ ■ - , s a ■ • ■ , s n }, H^V and (<, <) -> (>, >) in ©. 
One can also flip spin in other chains to get a new degenerate excitation. 



FIG. 1: (Color on line) The ground state ( taking \G^) as an example), and the low-lying excitations 
in a set of decoupled Ising spin chains which form a square lattice. From left to right and up to 
down, they are |Gf), \Hi a ), |^i a }, |Wp), |Wp P ,}, \W Ppl ) and |Wp/). Up- and down-arrows 

label the fermion with spin-up and spin down. Empty circle is unoccupied site and up-down arrow 
is double occupied. The white plaquette P" has (Gp„,Gp„) = (1,1), the yellow has (—1,1), the 
red has (1, —1), and the grey has (—1, —1). 

Fusion rules Note that O 2 = I for O = H, V, W, W d , W h . The fusion rules of these 
excitations are given in Table 1. The fusion rules of the closed subset {I ,V ia ,Wp,Wp P ,} 
(or {/, H ia , Wp, Wpp/} ) are exactly the same as the fusion rules in Kitaev toric code model 



if we identify 2\, Wp, and W PP , as ip,e and m (3, 9]. The subset {/, T ia , W Ppl , Wp P ,} 
has similar fusion rules but JF is bosonic. {/, 1~Ci a , T> ia , J-" ia } is also a closed subset and has 
similar fusion rules but with two fermions (H,T>) and one boson (J 7 ). 
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Table 1: Fusion rules of excitations. 

Discrete gauge symmetry The conserved quantities are simply Q s i = 2n s ^ — 1, (in fact, they 
are n s ^ ), which have eigen values ±1 at each site. This generates Z 2 x Z 2 gauge symmetry. 
For a plaquette P, one can label the plaquette by G S p = (2n iaySa — l)(2n JojSo — 1) for a pair of 
nearest neighbors (i a ,j a ) in the a-th chain, which also have eigen values ±1. Obviously, the 
ground state is of all Gp = 1. The excitations are also of all Gp„ = 1 except for the right 
plaquette P and left plaquette P' of i a where (Gp , G S p, Gp,, G%) = (-1,1,-1,1) for H ia , 
(1, -1, 1, -1) for V la , (-1, -1,-1, -1) for F ia , (1, 1, -1, -1) for W P , (±1, T l, T l, ±1) for 
Wpp, and (=f1, ±1, ±1, q=l) for Wpp,. To distinguish latter two, one uses (Q Sa ,i a ,Qs a ,ia) = 
(—1,-1) for Wpp, and (1,1) for Wpp,. In this sense, these string-like excitations are also 
called Z 2 x Z 2 vortices. 

Boundary conditions A string-like excitation may be thought as a spin flipping domain wall, 
a topological defect (See Fig. 1, late four). The open boundary condition allows odd number 
of domain walls while the periodic boundary condition restricts the wall number to even in 
each chain. If we curl up the plain lattice to a torus, the topological order is completely 
trivial and the system is identical to a one-dimensional one. In the open boundary condition, 
two walls as ends of a string, e.g, • • • ft II • • • and • • • |||t • • • , may locate at different 
chains so that one can circle around another. Therefore, their (mutual) exchange statistics 
can be well defined. 

Statistics in a given subset We now study the statistics of the non-local excitations within 
a subset. The local excitations H,V are the fermionic while T is bosonic. For a string-like 
excitation, the site of its domain wall, i a , may be used to label its 'position', e.g., Wp = Wj , 
etc. W is bosonic because it is a string of T . W d,h themselves are bosonic, e.g., 

Wfw£ = (W1XW1ZW2) = (H{D{H 2 )(Hx) = W 2 h Wf. 
Since W and W d,h are distinguishable, the exchange between them does not make sense. 
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However, because WW h,d fuses to a fermion while themselves are bosons, when W cir- 
cles around W h,d or vice verse, a minus sign is acquired. In general, this fact can be 
oroved by applying the consistent conditions, i.e., the pentagon and hexagon equations 



181 ] . For the present case, one can take Kitaev's graphical proof [9j]. For simplicity, denote 



{I,H,W,W d }( or {I,V,W,W h }) = {I,ip,e,rn}. Since e 2 



m 



I, we can create two 



pairs of e and two pairs of m from a given ground state at an initial time T Q (See Fig. [2]). 
These excitations move along the pathes as shown in Fig. [2j At certain time T±, they fuse 
to four fermions ip. As time flies, while blue and black fermions stay alone, green and red 
fermions exchange their positions. Finally, at T 2 , the fermions split to m and e pairs which, 
at the end (Tf), fuse back to the ground state. As two fermions exchange, this process con- 
tributes a minus sign = — 1 to the ground state comparing to a process without fermion 
exchange. Now, examine this process in non-local excitation exchanges. As shown in Fig. 
[31 this fermion exchange is corresponding to four exchanges: R em , R ee , Rmm and R me - That 
is. 



Rme Ree Rmm Rem 



Rme R< 



me ± ^em 



R 



-1, 



(4) 



since R ee = R mm = 1 as e and m themselves are bosonic. The minus sign when e and m 
doubly exchange, or equivalently, e encircles m, R me R em — — 1, proves the mutual statistics 
between e and m is semionic. Note that for the subset {/, W d , W h }, since T = W d yV h 

is boSOniC, Ry^dy^hRy^ihy^d (= -Rjc-jc- = 1) is trivial. 




nr. - m eve 



FIG. 2: (Color on line) Two pairs of m (dashed lines) and two pairs of e (solid lines) created from 
the ground state at To; they fuse into four fermions ip (colored ellipses) at Ti; then they split back 
to four pairs at T 2 and annihilate to the ground state at Tf. The green and red fermions exchange, 
which differ a minus sign from none of exchange. The arrow indicates the time direction. 



Peeling off a subset from degeneracy So far, we only say string-like excitations obey mutual 
semionic statistics but not call them anyons or semions. The reason for this is there are 
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FIG. 3: (Color on line) Divide the fermion exchange into exchanges between the same non- 
local excitations, R ee and i? mm , and double exchange e and m, R me R em , which is equal to moving 
e (or m) around m (or e). 

many degenerate non-local states which are not in the same subset. For example, nipping any 
chain's spin for W-p results in a degenerate excitation with W-p but already out of the subset 
of W-p. In this sense, these excitations can not be identified as individual quasiparticles. 

To peel off a designed subset, we need to set barriers between individual degenerate 
ground states without changing their energies. To control the electron spin of each individual 
chain is not easy. However, it becomes possible in a cold atom (molecule) system because the 
'spin' we are studying actually labels the different hyperfme states of atoms in the cold atom 
context . Once an atom is in a given hyperfme state, local fluctuations from the environment 
can not switch it to others. Therefore, we can peel off a given string-like excitation from 
others by preparing the ground state. For example, we can apply a magnetic field so that 
only the atoms with a given hyperfme state are loaded into the lattice and then turn off 
the magnetic field after the system is stable at the ground state. A global ferromagnetic 
ground state \G^) = IG^^...^}) is prepared. All excitations in (jHJ) then can be prepared 
by creating, annihilating the fermions or changing fermions from | to 1 hyperfme states by 
means of recently developed stimulated Raman spectroscopy or photoemission spectroscopy 
technique [3] . Here annihilating and creating a fermion do not means removing fermions 
from or reloading them into lattice sites. It can be turned into other hyperfme states which 
are almost not coupled to the 'spin' j and j hyperfme states or reverse. The non-local 
excitations may be controllably prepared. We may merely prepare excitations in a given 
subset and they are barricaded from their degenerate states. The non-local excitations 
obeying mutual simonic statistics in this subset are now identified as mutual semions 



Recently, a realization of this model in superconducting circuits has been proposed [20 ]. 
It is also a possible way to peel off a semionic subset. 
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Cold fermions with dipole-dipole interaction For cold fermions in optical lattice, off-site inter- 
action between the cold atoms can be induced by their diploe-dipole interaction. Recently, 



the degenerate Fermi gases of rare-earth atoms of Ytterbium (Yb) have been obtained [16 ]. 
They are possible candidate to be a practical system of our model because the fermionic 
isotopes 171 Yb and 173 Yb are stable in nature and their metastable state 3 P2 has a large 
magnetic dipole moment 3/ig. Deeply bound cold fermionic heteronuclear molecules have 
much larger dipole moment, e.g., the electric dipole moment of 40 K 87 Rb in its absolute bound 
ground state is 0.3 ea# [17J. Load the fermions in an optical lattice and polarize all dipoles 
along the horizontal diagonal of squares by using an external field. The interaction potential 
is Vrf(r, r') = d 2 1 ~ 3 ^,° ij2 6 where is the angle between R = r — r' and d (the diploe moment 
of an atom). The interactions along the diagonal become attractive. The repulsive interac- 
tion is restricted in the region with > O c for cos 2 C = 1/3. The interacting Hamiltonian 
can be written as 

{ij) hd ,s ((ij))hd,s ij,0<e<B c ,s ij,B>Bc,s 

where ((ij))hd denotes the sum along the horizontal diagonal other than the nearest neigh- 
bors. It is very easy to stabilize the ground state because one may increase the distance 
between the horizontal chains or adjust the optical lattice potentials so that the inter-chain 
couplings become weak. 

Strictly speaking, the anyons emerging from these dipolar particle systems are logarithmic 
confinement in thermodynamic limit, i.e., E pSuiT — E g oc InL as L — > oo for L the distance 
between a pair of the non-local excitations. This weak divergence in a practical optical 
lattice may be abided, e.g., if the short range model we proposed has -E pair — E g ~ 1, this 
logarithmic excitation energy is 4.5 for L = 50. Even L = 1000, this energy only increases 
about 15 times. 

Conclusions In conclusions, we have proved that it is possible to find a non-trivial mutual 
anyonic statistics in a fermionic system with conventional two-body interaction in open 
boundary condition. How to peel up a semionic quasiparticle from many degenerate states 
and manipulate them were discussed. 
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